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ABSTRACT 


The  steady-state  analysis  of  open  or  closed  networks  consisting 
of  n  nodes  interconnected  by  links  over  which  mass  flows  is  the 
subject  of  this  report.  Formulas  are  derived  which  yield  the 
equilibrium  mass  balance  at  each  node  in  terms  of  a  set  of  speci¬ 
fied  parameters  which  characterize  both  nodes  and  links.  The 
solution  is  then  generalized  to  the  case  of  a  k-priority  network  of 
interconnected  nodes,  either  open  or  closed,  k-priority  systems 
are  solved  recursively  and  formulas  are  given  which  determine  delays 
at  each  node  for  a  fixed  priority  mass  element  in  terms  of  known 


parameters . 
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INTRODUCTION 

A  dynamic  model  of  a  materials  balance  consisting  of  a  system 
of  simultaneous  equations  has  three  possible  modes  of  behavior 
as  time  increases.  There  may  be  unlimited  increase  or  decrease 
of  a  variable;  there  may  be  attained  constant  levels  of  the 
variables  not  depending  upon  time;  or  there  may  be  oscillatory 
solutions.  Thus  when  the  derivatives  of  the  dependent  variables 
are  set  equal  to  zero,  the  roots  of  the  resulting  system  of  equa¬ 
tions  identify  the  combinations  of  values  of  the  dependent 
variables  within  which  the  steady-state  combinations  are  a 
subset.  By  steady  state  one  means  either  the  case  in  which  the 
dependent  variables  attain  constant  values  independent  of  time 
or  the  case  in  which  they  yield  fixed  temporal  averages.  The 
linear  steady-state  network  model  of  an  open  system  referred  to 
here  can  best  be  visualized  in  terms  of  a  network  flow  diagram 
containing  several  Kodcc  or  points  of  accumulation  of  "species,’' 
i.e.,  substances  of  some  kind  connected  by  linkages  over  which 
material  transfer  occurs  that  is  governed  by  physical  laws  or  man 
made  control  rules  or  both.  A  basic  governing  principle  in  the 
construction  of  the  model  which  is  true  for  deterministic  cases 
and  for  a  wide  variety  of  cases  in  which  the  flows  are  governed 
by  stochastic  rules  [1 J  is  the  following: 


1 


total  means  input  rate  to  cell  i  = 

Q.  volume  of  cell  i 
1  _  _ _  s 

mean  delay  time  in  cell  i 

=  total  mean  output  rate  from  cell  i.  (1) 

s  identify  those 

inputs  and  outputs  of  substances  which  have  the  effect  of  forcing 
exchanges  among  cells  in  such  a  way  that  a  set  of  constant  flows 
across  nodes  results.  In  a  steady-state  condition  the  sum£e_^ 
necessarily  equals  the  sumEd^.  One  anticipates  that  this  condi¬ 
tion  can  be  satisfied  in  which  different  combinations  of  mean  flow 
rates  across  nodes  are  possible  and  in  fact  it  is  shown  below  that 
one  can  set  independently  the  values  of  mean  flow  rate  across  a 
fixed  subset  of  nodes  while  satisfying  the  side  condition =  £d^« 

The  solution  one  obtains  is  in  terms  of  the  constant  ratios 


A  set  of  external  sources 


and  sink 


which  are  the  limiting  values  obtained  by  averaging  the  input  rate 
to  node  i  over  time  from  0  to  +  «?,  Knowing  the  ratios  r^  and 
cither  the  mean  residence  time  D  or  cell  volume  the  remaining 
quant itiv  can  be  determined.  The  model  can  be  used  for  identifying 
feasible  sets  of  external  demands  upon  the  network  as  well  as 
obtaining  specific  solutions.  In  application  the  nodes  might 
represent  centers  in  which  traffic  of  some  kind  arrives  and  departs 
after  the  occurrence  of  certain  events  at  the  node.  The  sinks 


might  also  represent  demand  points  for  a  commodity  such  as  water 
over  a  geographic  area. 


EQUATIONS  Ql:  MASS  BALANCH 


An  open  system  differs  from  a  closed  system  in  that  there  are 
exogeneous  sources  which  feed  mass  into  the  system  at  one  or  more 
nodes  at  specified  rates.  In  addition  there  are  exogeneous  sinks 
into  which  mass  flows  from  certain  nodes  at  specified  rates.  If 
the  system  is  in  a  steady-state  condition,  the  mass  input  rate  must 
equal  the  mass  output  rate.  A  closed  system  is  obtained  if  this 
value  is  identically  zero.  A  specification  of  the  exogeneous 
inputs  e^  and  outputs  does  not  determine  the  flow  rates 
into  nodes  from  strictly  internal  deputations  within  the  system. 

It  does,  however,  restrict  the  set  of  possible  combinations  of 

and  hence  r^'s  that  will  satisfy  the  mass  balance  equations. 

By  definition,  one  has  in  steady-state  that 


D. 


1 


r . 

i 


+  e. 
i 


=  u.  +  d. 

i  i 


r . 

i 


where 


d.  =  mass  flow  rate  to  an  exogeneous 

sink  from  node  i  (i  =  1 , .  . .  ,n)  . 

e.  =  mass  flow  rate  into  node  i  from 

an  exogeneous  source,  (i  *  l,..,n) 

u.  =  non-negative  mass  flow  rate  out  of 
node  i  to  other  nodes  in  the 
network,  (i  =  1 , . . . .n) 

Although  the  average  through  puts  \r?l  must  exist  if  the  s vs ten 
has  a  steady  state^thoy  do  not  necessarily  represent  parameters 
exp  licit  i  i>  1  lie  f orrcu  1  a t  i  on  o  f  t  ho  t  i  me  d  eper.d  out  phy  s  i  ca  1  ayna; . i  o  s 
of  the  materials  balance  equations.  They  may  be  used,  to  construct 
steady-state  equations  of  balance  •  u  d  j  f  fert  nt  i  al  eqwu- 

tions  governing  tho  temporal  dynamics  of  the  system.  In  this  case, 


one  may  specify  independently  certain  direciiorial  coefficients 
which  are  needed  in  order  to  solve  for  the  which  yields 

rate  information  (not  accumulations).  Let 

P  =  [p4j]  (4) 

where 


i) 


=  fraction  of  mass  leaving  node 
i  that  arrives  at  node  j . 

(i  ,  j  =  1 , .  . .  * 


In  some  applications  there  may  be  limits  upon  the  rates  at  which 
mass  can  enter  or  leave  or  both  enter  and  leave  a  node.  Let 

a.  =  maximum  allowable  arrival  rate  of  mass  (5) 

1  into  node  i  from  combined  sources. 

(i  —  1 , . . ,  >n) 

a_  =  minimum  allowable  arrival  rate  of  mass  (6) 

into  node  i  from  combined  sources. 

(i  =  1 » » •  > n) 

b.  =  maximum  allowable  departure  rate  of  (7) 

mass  from  node  i . 

(i  -  1  >  •  •  • » n) 


b^.  =  minimum  allowable  departure  rate  of 
mass  from  node*  i  . 


Note  that  a.  =11.  +  e.  =  a.  and  b.  =  u.  +  d.  =  b,  (i  - 

-~i  i  l  1  —i  x  x  i  v  >  •  J 

For  the  i-th  node,  therefore,  one  has  under  steady-state  condition 
that 


or,  since 


M 

£ 

j-1 


u .  p  . .  +  e. 
i  3 1  i 


u .  +  d . 
(1=1,. ...n) 


(Sj 
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Equations  (9)  can  be  written  as  a  single  n  xn  matrix  equation 
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where 


Q  =  total  mass  conserved  in  system. 

In  applications  involving  materials  handling  systems  or 

certain  traffic  systems  the  p^ 1  s  may  be  fixed  in  advance  or 

it  may  be  clear  as  to  which  combinations  of  P—  's  arc  permissible. 

In  other  applications  where  one  is  dealing  with  large  volumes 

of  a  mass  such  as  lake  water  or  liquid  moving  through  a  porous 

medium  the  p..'s  mav  be  unknown  but  could  be  determined  as  part 
ri  1  "  1 

2 

of  the  solution  in  which  case  one  is  dealing  with  n~+n  unknowns. 
The  number  of  equations  in  this  case  depends  upon  boundary 
conditions  governing  the  and  other  physical  considerations. 


This  paper  deals  on 1 v  with  the  case  in  which  the  p. . ? 
*  4  '  *1.1 

known  or  can  be  specified  paranctically. 


s  are 


6 


rank  or  Tin:  p- matrix 


It  is  easy  to  show  that  the  n-th  row  of  the  P  matrix  in 

equation  CIO)  is  a  linear  combination  of  the  first  (n-I)  rows 

the  result  being  that  its  rank  is  at  most  n-1.  Thus,  at  least 

one  rate  ^  (for  some  k)  must  be  specified  independently  implv- 
Dk  .  . 

ing  that  many  different  internal  circulation  rates  are  possible 

while  the  external  inputs  ^  and  outputs  remain  fixed. 

However,  ^k  =  r.  may  be  limited  by  restrictions  (5) ,  (6),  and 
Dk 

(7),  i  .  e  . , 

^  Qk  o  -  - 

max  [o,  V  h]  =  rk  =  Uk  +  dk  "  =  min  lV  bkJ 


and  n 


(k  =  1 ,  .  .  ,n) 


Consider  the  three  node  network  specified  below.  Assume  an 

exogeneous  source  at  node  1  with  =  100  and  an  exogeneous 

sink  at  node  3  with  d„  =  100.  All  other  e.  and  d.  are  assumed 

j  it 

to  be  zero. 


Upon  solving  (13)  one  finds 


Q.  -  22  +  .37  Q, 

Q2  =  -46  +  .60  Q5  (14) 

In  other  words  the  system  has  one  degree  of  freedom  but  Q„  is 
constrained  by  inequality  (12).  Suppose  further  that 
a. ^  =  tr,  =  0  and  +  «  so  that  (12)  reduces  to 

o  =  u  +  100  =  ^3 
3  3 

o  %  u 

3 

where  upon 

Q.  =  300. 

If  =  300  then  «  0  meaning  that  node  5  acts  as  a  sink  rela¬ 
tive  to  internal  circulations  within  the  svstem.  If  one  sets  Q, 

>  300  then  iu  >  Q  and  node  3  no  longer  behaves  as  a  sink.  It 
is  also  clear  that  increasing  Q~  also  increases  the  internal  flow 
rate  between  other  nodes. 

It  is  apparent  that  provided  the  P  matrix  cannot  be  bloek- 
di agonal iced  to  the  form 


(15) 


meaning  that  all  nodes  "communicate the  rank  of  P  is  n-2.  If, 
in  contrast,  one  can  subdivide  P  as  shown  in  (15)  one  is  dealing 
with  independent  subnetworks  which  may  be  analyzed  separately. 

The  rank  of  P  in  this  case  is  where  f  is  the  number  of  blocks 
on  the  diagonal . 

MEAN  DELAY  IN  THE  SYSTEM 

The  mean  delay  (mean  residence  time)  at  a  node  is  related 

to  mean  flow  across  the  node  and  node  volume  according  to 

equation  (3).  The  values  of  any  two  of  the  quantities  u^, 

(L  ,  or  D.  determines  the  third  which  holds  for  deterministic 
i  i 

systems  and  for  many  stochastic  systems  [1,  2].  For  present 
purposes  it  is  assumed  that  mean  delays  at  nodes  are  known  or  can 
be  computed  once  the  u^fs  are  known.  The  mean  delay  of  an 
element  of  mass  in  the  system  can  now  be  computed  when  the  rule 
governing  the  transfer  of  mass  among  nodes  is  known.  Little 
more  can  be  said  unless  one  assumes  that  mass  elements  move 
within  the  system  at  random  so  that  at  a  particular  node  each 
element  has  an  equal  chance  of  being  selected  to  move  in  a 
particular  "direction."  In  this  case  the  theory  of  first 
passage  times  for  Markov  chains  applies  and  the  mean  delay 
of  a  unit  in  the  system  can  again  be  computed  in  a  straight 
f  o  rwar  J  m  a:  :ier .  The  c  omput  at  i  er.s  a  re  re  v  i  ewe  J  b  c  1  cw . 


M. 

1 J 


mean  dr  lav  encountered  by  a  mass 
elenr.it  moving  from  node  i  to 
node  j  .  (i  ,j  =  1  i  4}  ) 


Let 


The  value  of  M.  .  can  be  determined  in  the  following  wav  which 

is  a  straight  forward  application  of  the  theory  of  first  passage 

T 

times  lor  Markov  chains.  Denote  by  P  the  transpose  of  the  P 
matrix.  Let  j  denote  the  node  to  which  a  mass  element  is  to  be 
transported.  Redefine  the  p^'s  corresponding  to  node  j  as  (1, 

0,  - ,  0)  so  that  once  a  mass  element  enters  node  j  it  remains 

T 

there.  Reposition  this  vector  as  the  first  row  in  P  so  that  the 
T  — ' r 

revised  form  of  P  ,  call  it  P  is 


— T 
P 


0  —  0 


P  P  - P 

21  22  1 2N 


P  .  P  .  ---  P 
nl  n2  nn 


(17) 


The  only  requirement  on  the  order  of  the  remaining  rows 

and  columns  is  that  their  order  agree.  That  is,  row  order  must 

equal  column  order.  (In  many  cases  =  0).  Denote  by  S  the 

— T 

submatrix  obtained  from  P  by  deleting  the  first  row  and  first 

col unn.  Let  [in.  .]  be  the  matrix  defined  bv 
13 

In,..]  =  [I  -  S]"1 


Then 


M.  .  =  £-  m.  n. 

1 1  il  9. 


(18) 


Continuing  with,  the  previous  example,  and  assuming 

that  set  i  -  3  so  that 

ii 


10 


and 


Then 


and 


r 


T 

p  =  1 
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1/2  0 

1/2  1/2 


2 

0 

1/2 


S  = 


0  1/2 

1/2  0 


I  -  S  » 


[I  -  SJ 


-] 


The  mean  delay  encountered  In 


node  5  fro, 


r  a  rnass'  «lw.ent  in  reach in- 


no~c  1  is  therefore 


(20) 


(21) 
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MULTIPLE  PRIORITY  SYSTEMS 

Multiple  priority  systems  occur  when  mass  elements  are  tagged 
in  a  way  which  identifies  them  by  the  mean  delay  they  encountered 
at  a  node.  Whenever  one  "species"  takes  precedence  over  another 
at  a  given  node  in  such  a  way  that  it  encounters  a  shorter  time 
delays  it  is  said  to  have  a  higher  priority  than  the  other  species. 
The  only  complicating  factor  entering  into  the  analysis  of  the 
multiple  priority  system  is  the  computation  of  the  mean  delr*s  for 
the  lower  priority  elements  since  they  now  depend  upon  the  mean 
delays  of  the  higher  priority  elements. 


B  A!.  AN  Cl;  IX  A  2 -PRIOR  TTY  SYSTEM 

It  is  assumed  that  the  network  contains  two  "species"  of  mass 
which  flow  through  nodes  and  over  links  in  the  manner  described 
above.  The  difference  here  is  that  one  species  takes  precedent 
over  the  other  in  flowing  through  any  node.  The  "low  priority" 
specie^  flows  through  a  node  only  when  there  is  no  "high  priority" 
mass  at:  that  node.  Pcnoting  the  high  and  low  priorities  by  1  and 
J,  respectively,  one  defines  the  following  symbols. 


1 


Let 


L>  .  =  m~an  delay  at  node  i  of*  priority  1  mass  el  ement s 

li  J  ' 

(i  =  1,. . .  ,n)  . 


^2i  =  mean  delay  at  node  i  of  priority  1  mass  elements 


(i  = 


Q]i  ;  quantity  of  priority  1  mass  present  at  node  i 


under  equilibrium  conditions 
(i  =  1 , . . .  ,n) 


-  quantity  of  priority  2  mass  present  at  node  i 
under  equilibrium  conditions 
(i  »  1 


Further,  let 


and  denote  the  average  output  rates  of  mass  of  priorities 
1  and  2,  respectively,  from  node  i  (i  = 


:  =  proportion  of  priority  1  mass  leaving  node  .i 
that  is  directed  to  node  i. 


p^i .  =  proportion  of  priority  2  mass  leaving  node  i 
that  is  directed  to  node  j, 


(l  -  l,  2;  i 


where 


b,  .  =  maxim:;;  allowable  departure  rate  of  k-th 
k  1 

priority  mass  from  i-th  node. 

The  coefficient  p,  .  alwavs  lies  between  0  and  1  and  is  called  the 
ki 

load  fzoioi*  for  mass  of  priority  k  and  the  i-th  node.  This  means 
the  priority  may  change  for  a  fixed  mass  element  at  different 
nodes . 

The  following  assumptions  are  made  concerning  the  behavior 
of  the  system  containing  mass  of  both  high  and  low  priority. 

(1)  The  analysis  with  respect  to  priority  ]  mass  remains 
unchanged  from  that  given  above  since  it  does  not 
recognize  the  existence  of  priority  2  mass. 

(2)  If  priority  1  mass  is  removed  from  the  system,  then 
priority  2  mass  is  the  only  priority  pi'esent  and  the 
mass  balance  analysis  is  the  same  as  that  given  above. 

In  this  case,  the  can  be  independently  specified 
sin:e  there  is,  in  fact,  a  single  priority  to  be 
analvzed.  In  such  a  circumstance  denote  D„.  by  the  symbol 

2i 

L>*.  where 

2i 

V* .  =  mean  dolav  of  yrioritv  2  mass  at  node  i 

2i  r 

when  all  priority  1  mass  is  absent  from 
the  system. 

(3)  When  both  priorities  of  mass  are  present  the  mass  balance 

cuuat ions  for  priority  2  follow  t  lie  so  of  priority  1  after 

the  del  a  vs  have  been  determined. 

2 1 


II 


(51) 


The  delay  of  a  priority  2  mass  element  at  a  node  is  the  sum 
of  1)  its  delay  L>*^  when  no  priority  1  elements  ar»  at  node  i ; 

2)  the  delay  D  ^  due  to  the  presence  of  priority  1  mass  elements; 

3)  an  additional  delay  generated  by  the  arrival  of  priority  1  mass 

elements  during  the  intervals  D*  and  I),,  which  is 
°  2i  li 


Jli 


<D2i  +  Uii>: 


an  additional  delay  generated  by  the  arrival  of  priority  1  mass 

elements  during  the  interval  defined  by  P, .  (D*  +  D,  . )  which  is 

J  li  2i  li 


PH  <D5i  *  “li’- 


4)  and  so  on. 


When  summed,  one  has  that 


D2i  ■£ _  «ii  (“li  *  "n) 


i=o 


=  (D*  +  D  . ) 


■2i  li  l-pn 


(3 


Note  that  all  quantities  cn  the  right  hand  side  of  (32)  are  known* 
p  j  ^  is  obtained  f ram  the  e :: an t i  on 


'li 


Q  • 
xIi  lJ 


li 


An  equation  m.  ■ !  t.  i  ;  non  solved  in  which  all  par..'  ct  l"*s 

and  variables  pertain  to  priori  :y-2  ”j;-,  l  ’pen  >dul  ion  e:.c  obtains 

rlie  ratios  rv  ~  ^2  i  . 

hT 
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can  be  computed.  The 


From  (32)  one  obtains  D  .  and  thus  Q?. 
essential  point  is  that  the  mean  delays  for  lower  priorities  cannot 
be  specified  independently  as  in  the  case  of  the  priority-1  elements. 

BALANCE  IN  A  K-PRIORITY  SYSTEM 


One  can  now  analyze  closed  or  open  networks  having  k  priorities 

of  mass  units.  Assuming  that  lower  priority  units  do  not  affect 

the  behavior  of  higher  priority  units,  one  solves  for  (Q  ,  # . ,Q  ) , 

CQ21,...,Q2n)  > - .  %!»•••>  Qj.n)  respectively .  The  solution 

of  the  mass  balance  equations  for  the  i-th  priority  requires  the 

previous  (i-1)  solutions.  Definitions  that  correspond  to  (26) 

through  (30)  are  made  for  each  priority.  The  only  real  complexity 

that  enters  is  the  determination  of  D..  which  is  the  delay  of  a 

J1 

mass  element  of  the  i-th  priority  at  node  i.  One  can  show  that 
L>..  is  given  by  the  formula 

ji 


D..  =  D*.  +  D.  ,  . 
Jt  31  3-1,3 


L>  .  +  D  . 
2i  lx 


OO  A 

+  £  (<r  .  ,  f  n*.  +  l).  P_.  +  D,  .  ]  = 

J'1*1  Ji  3-1.1  2i  lx  1 


1-TT. 


J-l.l 


fD*.  +  P.  ,  .  + 
1 J  J  - 1 , 1 


+  IV  +  I) 
2j 


li 


(34 


where 


f-i..-' 


ii  +  p;i  + 


’j-l.i 


(j  =  „,k)  (i  =  1 .  ,n) 
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